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We show how center vortices and Abelian monopolcs both appear as local gauge 
ambiguities in the Laplacian Center gauge. Numerical results, for 5(7(2) and 
51/(3), support the view that the string tension obtained in the center-projected 
theory matches the full string tension when the continuum limit is taken. 



1 Introduction 

There are several approaches addressing the problem of understanding the 
mechanism of color confinement in non-Abelian gauge theories. The most 
popular share the idea that only a subset of the degrees ofpireedom are rel- 
evant for confinement. In the Abelian projection approacltra, one takes into 
account the maximal Abehan subgroup ?7(1)^~^ of the gauge group SU{N) 
and monopoles are tJoe effective degrees of freedom under study. In the Center 
projection approaclfl'a, it is the center group Zm which is considered and center 
vortices are the effective degrees of freedom under study. 

These two schemes are commonly believed to give alternative descriptions 
of confinement. Many analytical and numerical studies have been and are 
being performed using these two approaches. They give clear evidence that 
both Abelian and center degrees of freedom play a relevant role. 

The reduction of the gauge symmetry and the selection of the effective 
degrees of freedom is usually carried out by a partial fixing of the initial SU{N) 
gauge freedom. We show that, in the Laplacian Center gauge, monopoles and 
center vortices are closely related in a unified description and that the latter 
are the effective degrees of freedom relevant for confinement. Indeed center 
vortices are related to a more reduced gauge symmetry than monopoles. A 
fundamental issue of this study is the use of the Laplacian Center gauge, which 
is free from the problem of the lattice Gribov copies affecting the widely used 
Maximal Center gauge. 
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2 The center degrees of freedom 



Let us consider an SU (2) gauge field defined in a plane and expressed in 
terms of the radial and angular components (A^, A^). Suppose that the radial 
component is vanishing, A,. = 0, and the angular one is given by A^p = j^c^. 
This gauge field configuration describes a magnetic flux tube crossing the plane 
at the origin. If we consider a Wilson loop encircling the origin, we see that 
it has a non-trivial value e^'^'^^ — —1 with respect to the center group Z2 of 
SU{2). Conversely a Wilson loop non encircling the origin has a trivial value 
+1 with respect to Z2. This is an example of what a center vortex is. 

Consider now an SU{2) gauge theory on the lattice and decompose the 
link variable U^{x) as the product of two parts 

U^{x)^Z^{x)U'^{x) (1) 

where Zf^{x) lives in the center group Z2 and U'^{x) is the coset part belonging 
to SU{2)/Z2- For instance, this splitting can be carried out defining U'^[x) 
as having positive trace. If W{C) is a Wilson loop along the closed path C, 
making use of the decomposition (|^), we can write 



W{C) = a{C) W'{C) 



n ''^P) 
pes 



W'{C) (2) 



W'{C) and cr(C) = Hpes '^^P) ^"^^ respectively the Wilson loops evaluated 
with the coset links U'^{x) and with the center links Z^{x). IlpGi; 
product over all the plaquettes p belonging to a surface E supported by C; 
the value of cr(C) does not depend on the choice of S and, for simplicity, we 
can choose it as the planar surface bounded by C. If we fix a gauge where 
U'^ is smooth, then W{C) has a non-trivial value with respect to Z2 if cr(C) 
does. Moreover, considering (||) for a single plaquette, ct(p) = — 1 is a "signal" 
for a center vortex. In 4-dimensional space-time, center vortices form closed 
surfaces in the dual lattice. 



It is possible to give a handwaving and qualitative argument showing how 
the described center vortices can give rise to an area law behavior for the Wilson 
loops. According to the previous argument, we are interested in the plaquettes 
p having cr(p) — —1. Suppose that t is the probability for a plaquette to have 
a{p) = —1, consequently {l — t) is the probability to have a{p) = +1. Consider 
a time-slice of a lattice SU{2) gauge configuration, then center vortices form 
1-dimensional closed strings in the dual lattice. The following figure displays 
an example of such a time-slice 
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The drawn plaquettes have a(p) = — 1 and the broken hne is the closed 
string in the dual lattice representing the center vortex. Consider a Wilson 
loop W{C) which, in the figure, is represented by the continuous line: the 
plaquette in bold line belongs to the planar surface S spanned by C. Let L 
be the linear extension of C and ^ the average linear size of the center vortex 
strings. If ^ ^ _L, the value <7{p) of each plaquette p in S is independent of 
the others; then, with respect to the center degrees of freedom, one can write 

< W{C) >~< a{p) >^= e^ i°s(i-2*) (3) 

where A is the area of E. With this qualitative argument we have obtained 
on one hand that the center degrees of freedom can give rise to an area law 
behavior for the Wilson loop and, on the other hand, that the string tension 
is about twice the probability t for a plaquette to have a{p) = — 1. We have 
assumed that L: in order for this inequality to be satisfied for arbitrarily 
large Wilson loops, ^ must be divergent, that is center vortex strings must 
percolate through the lattice. Conversely, suppose that the center vortices do 
not percolate, so that ^ has a finite value, then for L sufficiently large, one 
must have ^ <^ L. The next figure displays such a case 




W(C) 



L 
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The values cr(p) of the plaquettes in E weU inside the boundary C are 
strongly correlated by pairs; moreover these pairs do not give a net contribution 
to cr(C). So only the plaquettes belonging to a ring of thickness ^ around the 
boundary of E can give a non-trivial contribution. According to this reasoning 
and, with respect to the center degrees of freedom, we can write 

< W{C) >^< a{p) >«^= e«^'°s(i-2t) 

where P is the perimeter of E. Thus, we have obtained a perimeter law be- 
havior for the Wilson loop. The conclusion from this handwaving argument 
is that the center degrees of freedom can give rise to an area law for the 
Wilson loop and that the area/perimeter behavior can be recast in terms of 
percolation/non-percolation of center vortices. 

3 Laplacian gauge fixing 

Many lattice-jSimulations - starting with the initial results by Greensite and 
collaboratorfl - have been performed in order to investigate the role of the 
center degrees of freedom in the non perturbative features of the non-Abelian 
gauge theories. In these studies the selection of the center degrees of freedom 
is carried out by the numerical partial gauge fixing of an ensemble of config- 
urations. The most widpJii used gauge is the Maximal Center gauge. In this 
gauge, it has been showiQ'Ll that the removal of the center vortices leads to the 
loss of the area law for the Wilson loop, to the restoration of chiral symmetry 
and to the disappearance of non-trivial topological features. The numerical 
implementation of the Maximal Center gauge fixing is performed by a local 
iterative procedure and, as there are niany local extrema, lattice Gribov copies 
are present. This is a serious problenU'Ll that can lead to a complete loss of 
meaningful information in the Z2 projected model. Thus, it is important to 
study the role of the center degrees of freedom in the confinement mechanism 
considering a smooth gauge not affected by this problem of the lattice Gribov 
copies. In the Laplacian gauge the reduction of the gauge degrees of freedom 
is carried out in an unambiguous way. This gauge was proposed by Vink and 
WiesoEl: they suggested to use the eigenvectors of the covariant Laplacian op- 
erator to fix the gauge. Since we are interested in reducing unambiguously 
the symmetry of the gauge group from SU{N) to its center Z^, it is useful 
to consider the Laplacian operator in the adjoint representation. In fact, as 
the adjoint representation is invariant under gauge transformations in Z^r, the 
adjoint Laplacian procedure fixes unambiguously the gauge up to the center 
symmetry. 
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Consider the 4-dimensional lattice SU{N) gauge theory. The adjoint co- 
variant Laplacian operator Ay'^{U) is given by 

- Af,{U) = ^ (2 6y,, 5-' - U;\x - A) 5y,.--, - Ul%^) 5y,.+ -^ (5) 

where a,& = 1, . . . , [N"^ — 1) are color indices and x,y are space-time lattice co- 
ordinates. The dotted U^{x) are the link variables in the adjoint representation 
and are related to the links U^{x) in the fundamental by 

t>f(2^) = ^Tr(Aa[/^(x)AbC/t(x)) (6) 

Ai, i = 1, . . . , (iV^ — 1) being the generators of SU{N) with the normalization 
Tr(AaAf,) = 25ab- If V is the volume of the lattice, A({/) is a [(iV^ - 1)^] x 
[(A^^ — real symmetric matrix which depends on the gauge field. The 
eigenvalues Hj of A are real and the eigenvector equation is 

Af^ijJWix)^^,,cp^^\y) (7) 

where <p'^^\ j = 1, . . . , [(iV^ — 1)V] are the (real) eigenvectors. So we can asso- 
ciate {N"^ — 1) -dimensional real vectors (j)^^^ {x) to every lattice site x. Consider 
now a gauge transformation on the links in the fundamental representation 
U'^{x) = ^{x)U ^^{x)^!^ {x + /i); the eigenvector equation becomes 

h'^'^\y)/^^^{U')(l^\x)cp^^\x)=^,,<j,^^Hy) (8) 

where VL is the gauge transformation in the adjoint representation. This re- 
lation shows that the eigenvalues are gauge invariant and the eigenvectors 
transform according to Cl°-'^{x)(l)jii\x) — (j)'^\x)' . This transformation law can 
be rewritten as follows 

Vt{x)<^^^\x)n\x) ^ <^^^\x)' (9) 

where we have defined the su{N) matrices - i.e. in the SU{N) algebra - 

= Eali^^aHx)Xa and = Eali'^aHxyXa. Gaugc trans- 

formations rotate the vectors (a;) in color space and we can fix the gauge by 
requiring a conventional arbitrary orientation for the 4'^^\x). To perform the 
reduction of the gauge symmetry from SU{N) to Zn we only need to fix the 
orientation of two eigenvectors of the Laplacian operator. As we are interested 
in the non-perturbative features and in fixing a smooth gauge, we consider the 
two lowest-lying eigenmodes 0*^^^ and c/)'^' . 
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The gauge fixing procedure can be split into two steps. In the first, one 
rotates ^^^\x) at every x so that ^^^\xy is diagonal. This leaves a resid- 
ual symmetry corresponding to gauge transformations in the Cartan subgroup 
U{1)^~^. To further reduce the gauge freedom we have to consider a second 
step where the second eigenvector (j)^"^^ is taken into account. The gauge trans- 
formation that has rotated to the Cartan subalgebra, maps to 
$(^)(a;)'. The remaining C/(l)^~^ symmetry can be fixed to by requiring 
that some conventionally chosen color components of the twice rotated ma- 
trix vanish. Now we describe explicitlv how to perform the presented 
two-step program for the SU (2) gauge theorjH, where, at each x, we consider 
the 3-dimensional real vectors and ^^^•'(a;). 

Step 1. We consider equation (||) for the first eigenvector and we con- 
ventionally define fl{x) to be the gauge transformation that rotates (j)'-^^(x) 
along the direction 3 in color space: i7(a;)0^^^ (a;) = cj)^^\xy oc (0,0,1). fl(x) 
is unambiguously defined up to gauge transformations V{x) in U{1). So we 
have reduced the gauge synmi£±ry to the Cartan subgroup of SU{2). This is 
the Laplacian Abelian gaugfiilO. 

Step 2. We apply the gauge transformation U{x) found in Step 1 to 
(jj^^^x): = fl{x)(j}^'^^x). 4>^'^\x)' is not invariant under the rotations 

V{x) £ U{1) that leave (j)^^\xy unchanged; then we can fix this symmetry by 
requiring, for example, that = V{x)(j)^^'> {xY lie in the 1-3 color plane 

in the positive direction. 

We have completely fixed the gauge symmetry in the adjoint representation 
and, as it is center-blind, we are left with the center^ymmetry The 
described two-step procedure for 5/7(2) can be extendectj to SU{N) with the 
same pattern of gauge symmetry fixing: SU{N) [/(l)^^^ in the first step 
and U{1)^~^ Zn in the second one. 

Two last, important remarks concern the accidental degeneracy of fii and 
fi2 in (0) and the arbitrariness in the eigenvectors If it happens that fii or 
fi2 is degenerate, the gauge fixing can not be carried out unambiguously and 
one has a (global) Gribov ambiguity. This case is really exceptional and never 
occurs in the numerical simulations. The second point is about the scale and 
sign arbitrariness of the eigenvectors 0^-'^ . Rescaling can not give rise to any 
ambiguity in the procedure while the freedom in the choice of the sign does. 
This global freedom can be eliminated with a conventional choice on 0*^^^ . 

4 Local gauge ambiguities 

The adjoint Laplacian gauge fixing procedure has local defects. Now we discuss 
how these defects can show up and how they can be identified with monopoles 
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and center vortices in SU{2). This discussion can be generalized to SU{N<^. 

Step 2 ill— defined: the second step of the Laplacian gauge fixing is not 
defined at the points x where (jy^^^'lx) / / 4>^'^^{x). In such a case also is 
invariant under rotations V{x) e U{^)- The condition 0*^^^(a;) // 0'-^-' (x) sets 
two constraints and so these points x - where the gauge symmetry is promoted 
from Z2 to C/(l) - form 2-dimensional surfaces in 4-dimensional space-time. 

Step 1 ill— defined: the gauge fixing procedure can not be even started at 
the points x where 4''^^\x) = (0,0,0). These defects constitute 1-dimensional 
strings in the 4-dimensional space-time since 3 constraints must be satisfied. 
At these points the symmetry is not fixed and the gauge freedom is SU{2). 

The 2-dimensional surfaces of Step 2 ill— defined can be identified with 
center vortices. Suppose that at a point xq it happens that (jP-\xQ) / / t/)^^-* (xq), 
then moving along a small loop around the singularity point xq, (/>^^^ describes 
a 27r rotation in color space. As a 27r phase in the adjoint representation of 
SU{2) corresponds to a tt phase in the fundamental one, it follows that a 
Wilson loop encircling Xq has a non trivial value with respect to the center Z2. 

The 1-dimensional strings of Step 1 ill— defined can be identified with 
monopole world-lines. In fact, in analogy with the Georgi-Glashow model, 
the points x where the Higgs field vanishes and the gauge symmetry can not 
be reduced from SU{2) to C/(l), correspond to the monopole world-lines. 

In the adjoint Laplacian gauge, monopoles and center vortices turn out to 
be closely related in a unified description. Consider the 2-dimensional surface 
S of the center vortices. At every x G 5, 0(^)(x) and (/)^^)(x) are parallel or 
anti-parallel. So S is divided in parts where 0(^)(x) is parallel to ^^^•'(x) and 
parts where it is anti-parallel. By continuity, these patches must be separated 
by 1-dimensional strings where (/)'^^^(x) = or (f>^^\x) — 0. Moreover, in the 
neighbourhood of a monopole, ^'•^■'(x) has a hedgehog-like shape and there will 
be a direction where 0'^-' // Thus monopole world-lines are embedded 

within the 2-dimensional surfaces of center vortices. 

5 Numerical results and their interpretation 

We have performed numerical simulations to investigate the role of the cen- 
ter degrees of freedom in the SU{2) and SU{3) lattice gauge theories. For 
SU{2) we have collected 1000 configurations at /3 = 2.3, 2.4, 2.5 on a 16'* 
lattice; for SU{3) we have generated 500 configurations on a 16'* lattice at 
/3 = 6.0. The following figure shows the measurement of the Creutz ratios 
X{R) = -ln{{W{R,R)){W{R-l,R-l))/{W{R,R-l)y) for 5t/(2) at /? = 2.4 
{W{R, T) is the R x T Wilson loop). Crosses refer to SU{2), circles to center 
projection after Laplacian gauge fixing and stars to the coset part. The contin- 
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uous band is the value in the hteratur(Bll3 for the SU (2) string tension at the 
considered set of parameters. The results show, on one hand, the flattening 
of the Creutz ratios in the Z2 sector and, on the other hand, the vanishing of 
the Creutz ratios computed with the coset links. We have obtained a similar 
behaviour for (3 = 2.3 and 2.5. In the case of SU{2>) also, the following figure 
shows the Creutz ratios in the sector after Laplacian gauge fixing. The 
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continuous band is the value in literatures for the 5[/(3) string tension at the 
chosen set of parameters. Also in this case, one can clearly see flattening to a 
non vanishing value for the Creutz ratios evaluated with center projected links. 

The good agreement with the values in the literature for the string tension 
in SU{2) and SU{3) should not be over-estimated. Numerical simulations are 
performed at finite lattice spacing and lattice artifacts can give non-negligible 
effects in the center projected theory. Our conjecture is that even if, at finite 
lattice spacing, the flattening value of the Creutz ratios in the center sector 
changes with the particular lattice Laplacian used to fix the gauge, this depen- 
dence vanishes in the continuum limit. The observation of such a behaviour 
would be a robust confirmation of the relevance of the center degrees of freedom 
in the confinement mechanism. To investigate this issue, we have considered 
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three different lattice Laplacians differing by irrelevant operators: they have 
been built using smeared links in (|6|). Every set of 1000 configurations at 
P = 2.3, 2.4, 2.5, has been fixed in each one of the three gauges and the Creutz 
ratios have been measured after center projection. The table summarizes our 
results: Ri — y^(Ji/crsu{2) where cr,; is the string tension measured in the Z2 





13 = 2.3 


f3 = 2.4 


/3 = 2.5 


Ro 


0.813(23) 


0.860(20) 


0.978(18) 




0.592(12) 


0.720(11) 


0.804(12) 


R2 


0.547(8) 


0.653(7) 


0.739(11) 



sector; i = 0, 1, 2 is an index for the three lattice Laplacians and crgjj(^2) is the 
value in literature for the SU (2) string tension at the three values of (3. Thus, 
it is in the continuum limit (/3 — > 00) that the string tension measured after 
center projection correctly reproduces the value of the full gauge theory. 
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